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The origin of observed extremely high energy cosmic rays remains an astrophys- 
ical enigma. We show that a single evaporating primordial black hole should 
produce 8.5 ■ 10 14 particles over a 10 20 eV threshold. This emission results from 
direct production of fundamental constituants and from hadronization of quarks 
and gluons. The induced flux on the Earth is studied as a function of the local 
density of exploding black holes and compared with experimental data. The dis- 
covery potential of future detectors is finally reviewed. 
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1 Introduction 



Small Primordial Black Holes (PBHs), with masses well below the self-gravitational 
collapse limit and possibly as low as the Planck Mass (Mpi pa 5.5 x 10 -5 g) may 
have formed in the primordial universe Jl| . Numerous processes, compatible with 
standard cosmological scenarios, can be put forward to explain their formation 
[0 . In particular, if they result from initial density perturbations (with an initial 
mass determined by the horizon mass at this epoch), the mass spectrum can 
be analytically determined [§] following the natural hypothesis of scale-invariant 
Gaussian fluctuations: 

(a - 2)M~ a M:- a lQ PBHPcrit (1) 



dMd y V / evap 

where M evap pa 10 15 g is the mass of a PBH evaporating nowadays, Vt PBH is the 
current density of PBHs in units of critical density p crit and a = (1 + + 
7) + 1, 7 = p/p being the pressure to density ratio. 

This study is dedicated to the final-stage emission of PBHs to investigate 
if they can be considered as candidates for extremely high energy cosmic rays 
(EHECR), beyond 100 EeV (10 20 eV). Observational data U show that the cos- 
mic ray flux seems to be curiously unaffected by the expected Greisen-Zatsepin- 
Kuz'min (GZK) cutoff (due to interaction with the 2.7K cosmological background 
above photoproduction threshold). The integrated emission of PBHs is estimated 
in the following sections for a volume of universe where predicted effects of this 
interaction are weak (i.e. for a radius close to the attenuation length pa 50 Mpc). 
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2 Individual emissions 

Hawking showed || that black holes can radiate particles in a process qualita- 
tively equivalent to e + e~ pairs production in a strong electric field. When the 
hole temperature becomes greater that the quantum chromodymanics confinment 
scale (T > Aqcd), i-e. some hundreds of MeV, emitted particles are fundamen- 
tal constituents rather than composite hadrons ||. The EHECR production by 
PBHs (which are particularly affected by evaporation effects because of their low 
mass) has to be understood in such an approach. The emission spectrum for 
particles of energy Q per unit of time t is: 

where contributions of angular velocity and electric potential have been neglected 
since the black hole discharges and finishes its rotation much faster than it evap- 
orates @ ||. k is the surface gravity, s is the spin of the emitted species and 
T s is the absorption probability. In the general function of Q, the 

particle mass m, the hole mass M, and the number of degrees of freedom of the 
species. Its value can only be computed by numerical approximations based on 
expansion in spherical harmonics of the scattering matrix || . In the optical limit 
(i.e. Q — ► oo) , which is totally justified for energies considered in this work, 

. 27Q 2 



3 6An 2 (kT) 2 
where T is the "temperature" defined by 



(3) 



In such a description, the black hole behaviour mimics a black body whose tem- 
perature increases when the mass decreases until it reaches the Planck limit where 



this theoretical description becomes unadapted. The time-evolution of the system 
depends on the emitted constituents' degrees of freedom and is therefore based 
on the choice of a particle physics model. It is likely that new particles, absent 
from the standard model, are emitted when the black hole temperature becomes 
extremely high, but the general behaviour remains unchanged: all the emission 
above 100 EeV is nearly instantaneous. The mass loss rate of a PBH is [|Kj [|ll 

dM _ (7.8d s= i /2 + 3.1rf s =i) ■ 10 24 gV 1 



dt M 2 



(5) 



where d is the mass- dependant number of degrees of freedom for the emitted 
particles of spin s. In the standard model, dM/dt pa —7.9 x 10 26 /M 2 above the 
top quark production threshold. It leads to 

1 h 3 c 6 d(kT) 

1 ~ (7.84=i/2 + 3.14=i) ' (4tt)6G 3 ' JkJf U 

or 

where t* = t/ls and kT* = fcT/lEeV. Since it has been checked that only 
particles emitted at kT > 5 EeV will contribute (within a few percent) to the 
flux of cosmic rays with energies beyond 100 EeV, the characteristic production 
time is At < 4 x 10~ 18 s. As a comparison, the total evaporation time for a 10 15 g 
black hole is of the order of the age of the universe. 

3 Extremely high energy emission 

Taking into account formula (6) relating the temperature to the mass, the pre- 



vious emission spectrum can be rewritten in its integral form |T3 per particle 
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species above a threshold E t h- 

1 27h 2 c 9 r kT pi 1 r°° Q 2 d(kT)dQ 



N 



(7.84=i/2 + 3.14=i) 8 6 7r 8 G 3 Jm {kTf J Eth eQ/( fcT ) - (-l) 2s v/ 
where Tj and Tpi are the initial and Planck temperatures. It can be numerically 
expressed as: 

N = 1 56 ■ 10" ^ r ^ (9) 

At« Wt) e* - (-1)2- W 

where x = Q/(kT). Fig [1] shows that the Planck cutoff is effective for energies 

well beyond those of interest here. 

After their production, emitted quark and gluons fragment and produce a 

subsequent number of hadrons. Monte Carlo simulation codes tuned to reproduce 

experimental data obtained on colliders cannot be used because the energies 

considered here are several orders of magnitude greater than those available today. 

The multiplicity of charged hadrons produced in a jet of energy Q is therefore 

estimated by means of the leading log QCD computation [ [l"3|j : 



n h (Q) « 3 x io-W^Q/a) + 2 (10) 



To get the resulting hadron spectrum, Hill |14| derived the following distribution: 



dn h , o v . nil (1 - z^ 2 



-± « lO-WW** 1 . x ^—=4- (11) 
dz z ^ lr ^z 

The number of emitted hadrons above the threshold, by a PBH of temperature 
T, can then be written as: 

rkTpi„ d(kT^) f°° x 2 dx f 1 dn h 



N h = 1.56 • 10 ib / ,) :i / , — / — dz (12) 

JkTi, \kT*Y Jm&iikT) e x — (— l) 2s JE th /(xkT) dz 

per particle species of mass m. The numerical computation has been compared 
to what is given by the empirical function [IT] drih/dz = (15/16) x z~ z l 2 {\ — z) 2 , 



leading to a multiplicity which can be easily calculated analytically. Results are in 
agreement within an error of 12% which is certainly not the dominant uncertainty 
in this evaluation. Figure |2] illustrates the general behavior of the total hadronic 
multiplicity fE th /Q ^JT^ Z above a given threshold. 

The total number of emitted particles above a detector threshold E t h = 
100 EeV can then be estimated by summing the direct flux (taking into ac- 
count all the standard model degrees of freedom) of fundamental stable particles 
and the fragmentated flux resulting from the previous computation for coloured 
objects. The numerical result is F(> 100 EeV) ~ 8.5 ■ 10 14 particles over the 
lifetime of a PBH. 

4 Resulting flux above 100 EeV 

The derivation of the exact resulting spectrum on the Earth is a complete study 
by itself, well beyond the scope of this work. It is straightforward to demonstrate 
that the integrated emission goes as E~ 2 , which seems quite difficult to conciliate 
with the cosmic-ray experimental data if energy- dependent confinement effects 
are ignored. The following section therefore aims at evaluating the orders of mag- 
nitude involved. 



To derive the resulting flux reaching the earth, PBHs have been considered as 
classical (non baryonic) cold dark matter clustered in galactic halos. The Milky 
Way mass distribution is therefore assumed to follow the simple law in spherical 
coordinates: 

^=^ir# (13) 
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where R is the distance between the considered PBHs and the Galactic Center, 
Pq is the local density of exploding PBHs, and R c is the core radius of the halo. 
The particle flux becomes: 



^jj-J = Po x F x J(R h ,R c ,Rq) 

s / galactic 



(14) 



where 



1 rn i-Rh 2 Rl + Rp 

J (R H , R c , R&) = — / R St 
07T Jo Jo Hi 







sincj) dR d<p 



(15) 



R\ + R 2 R%- 2RR Q cos<p + R? 
Rh being the total radius of the halo. Table |l] gives fluxes normalized to the 
average for extreme values of R c and Rh (for R & = 8 kpc): it shows a quite low 
dependance on the halo parameters. 



Relative fluxes 


R H =A0 kpc 


R H =100 kpc 


R H =200 kpc 


R c =2 kpc 


0.97 


1.01 


1.03 


.R c =4 kpc 


0.94 


0.99 


1.01 


R c =6 kpc 


0.96 


1.02 


1.05 



Table 1: Relative fluxes for different halo parameters 

The extragalactic contribution is computed by assuming a standard galaxy 
distribution pa ~ 0.01e ±0,4 /i 3 Mpc -1 [0 (with the Hubble parameter defined as 
Ho = h ■ 100 km.s _1 .Mpc _1 ). The resulting flux is: 



dAf 
dt 



K(F, R H , R c , Rq) x p Q x p G x R GZK 



(16) 



extragalactic 

where K(R H , R c , Rq) is the average emission of a single galaxy (obtained by the 
previous method) and Rqzk ~ 50 Mpc is the radius of a sphere "unaffected" by 
the GZK cutoff. On such distances, it is not necessary to redshift energies within 
the expected accuracy of a few percent. 
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Numerical results for average values of physical parameters show that the 
galactic contribution is nearly three orders of magnitude larger that the extra- 
galactic component, even assuming the highest galaxies number density and the 
upper Hubble parameter limits (h<l). As it only depends linearly on Rgzk, the 
accurate determination of this radius is also irrelevant. The total flux above 100 
EeV is then: 

(iS) = 3 - 8 ■ 10 " 23 >< f i 3 ) m-W" 1 (17) 

\ dt J pbh V 1 y ear p c / 



exp 



Experimental data on EHECR show an integrated flux of the order of (j^f) 
1Q-16 m ^ 2 s _1 sr _1 ||17|| . The required density of exploding PBHs near the earth 
to reproduce such a signal is then p « 2.6 ■ 10 6 year _1 pc~ 3 . 

5 Discussion 

Direct observational constraints on the local PBH explosion rate p & are quite 
difficult to obtain. A reliable search for short bursts of ultra high-energy gamma 
radiations from an arbirary direction have been performed using the CYGNUS 



air-shower array ||18|| . No strong 1 second burst was observed and the resulting 
upper limit, based on the exhaustive analysis of a very fine binning of the sky, 
is in the range p & < 0.9 ■ 10 6 year _1 pc -3 . Very similar results were derived by 



the Tibet [|T^] and the AIROBIC collaborations p0 |. TeV gamma-rays have also 
been used to search for short time-scale coincidence events, thanks to the imaging 
atmospheric Cherenkov technique developped by the Whipple collaboration. The 
very high-energy gamma-ray bursts detected are compatible with the expected 

8 



background, within ±1.7er. The resulting upper limit obtained with 5 years of 
data PHI , i.e. p < 3 • 10 6 year -1 pc~ 3 , is substantially better than the previous 



published results in the TeV range. All those limits are roughly compatible with 
the density required to generate the observed EHECR spectrum. 

At the opposite, low-energy (< 0.5 GeV) cosmic-ray antiprotons detected by 



a BESS 13-hours ballon flight have been used |L5] to put a much more severe 
upper limit of p© < 2 • 10~ 2 year _1 pc -3 , which could exclude PBHs as serious 
candidates for EHECR. This analysis is particularly promising since the authors 
have shown that the local PHB-antiproton flux can only be due to contributions 
from black holes that are very close to explosion, and exist within a few kpc away 
from the Solar system. Nevertheless, such data suffer from an important lack of 
statisics and from contamination effects due to interactions with the atmosphere. 
Future results from the AMS spectrometer on board the International Space 
Station will give a much more accurate antiproton spectrum in the 0.1-1 GeV 
range. Those data should allow a stringent upper limit (if not a positive detec- 
tion) on nearby exploding PBHs. 

An entirely different approch is to study the diffuse gamma-ray background 
spectrum. The emission from PBHs over the lifetime of the Universe is integrated 
so as to evaluate the resulting particles and radiations. This method |TT[ leads to 
Pe < 10year _1 pc -3 for clustered black holes. It should, anyway, be emphasized 
that such a study does not directly constrain p & . The resulting "Page-Hawking 



bound" |23| on VLpbh, derived to match the observed spectrum at 100 MeV, is 



converted into an upper limit on the initial number density of holes per loga- 
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rithmic mass interval Npbh at M = M evap under assumptions on the Hubble 
parameter, on the relative matter density (fijif), on the equation of state of the 
Universe at the formation epoch, and on the gaussian distribution of initial den- 
sity perturbations. This latter point is rather controvertial. The upper limit 
on po which can then be derived has to account for the large (possibly up to 
8 orders of magnitude) uncertainties associated with clustering. Recent reviews 
on the detection of PBHs captured around massive objects p4 show that, when 



the first astrophysical objects with masses of the order of the Jeans mass were 
forming, black holes haloes on the sub-galactic scale could have formed around 
old globular clusters, dark matter clusters or population III stars. This makes 
the use of 100 MeV gamma-rays a quite difficult way of ruling out an important 
local rate of PBH explosions, though future GLAST p7| data should change the 



situation by a dramatic improvement in statistics and resolution. 



Furthermore, the first results from the AUGER Observatory p8] will soon give 



high statistics samples of EHECR. With the PBH space distribution previously 
assumed, the resulting EHECR flux would be from 6.0 to 2.2 times higher in the 
Galactic Center direction than in the opposite direction, for an integrated obser- 
vation angle from 10 to 90 degrees. After five years of operation, the AUGER ob- 
servatory should collect up to 300 cosmic rays above 100 EeV. Such an anisotropy 
would be detectable if more than approximately 50% of them come from PBHs. 

Finally, some evidences for PBH signatures available nowadays should be 
noted. Studies of the BATSE IB and 1C catalogs have shown PSl that some 
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gamma-ray bursts (GRBs) were consistent with a PBH evaporation origin at the 
quark-gluon phase transition. Characteristics of selected events are in remarkable 
agreement with the "Fireball" PBH picture. The resulting (model dependent) 
limit is significantly lower than what is expected in the present work, and dis- 
favours a PBH origin for EHECR. Nevertheless, new analysis of EGRET data 



2q| gives some evidences for a gamma-ray halo "glow" due to PBH emission. 
Those first tentative detections are very promising for further investigations on 
the subject. 



From the theoretical point of view, it should also be emphasized that re- 
sults given in this paper are based on the standard particle physics model. The 
probable increase of degrees of freedom available when the black hole tempera- 
ture exceeds energies currently available on colliders would modify the estimated 
fluxes, making the final explosion much more violent. This could validate PBHs 
as a good source candidate for a fraction of the observed high energy cosmic rays. 

Acknowledgments. I would like to thank Cecile Renault for very helpful 
discussions. 
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Figure 1: Planck cutoff effect on the integrated primary spectrum 
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Figure 2: Hadron multiplicity as a function of the detection threshold and of the 
jet energy 
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